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Basic peculiarities of market price fluctuations are known to be well described by a recently 
developed random walk model in a temporally deforming quadric potential force whose center is 
given by a moving average of past price traces [Physica A 370, pp91-97, 2006]. By analyzing 
high-frequency financial time series of exceptional events such as bubbles and crashes, we confirm 
the appearance of nonlinear potential force in the markets. We show statistical significance of its 
existence by applying the information criterion. This new time series analysis is expected to be 
applied widely for detecting a non-stationary symptom in random phenomena. 
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I. INTRODUCTION 

Financial bubbles and crashes have been occurring oc- 
casionally in nearly every market causing social troubles 
of various magnitudes. In the ordinary market state, 
prices fluctuate fairly randomly and directional predic- 
tion is almost impossible; however, during bubbles and 
crashes, the market prices are known to move quite asym- 
metrically Viewing the movement of prices in a 
coarse-grained way reveals that the growth of price in 
a bubble period is approximated either by an exponen- 
tial function 0, [H , a double exponential functionjjL or a 
function having a finite-time sing ularity @, 1, 0, EM O • 
In any case, it is a reasonable assumption that there are 
some special mechanisms of bubbles that are quite differ- 
ent from the ordinary market state. It is very important 
to develop the teqnique to quantify the financial risk of 
large price changes 11 1 . 

In the year 1900, Bachclicr introduced the first random 
walk model as a model of market price fluctuations [l2| . 
The stochastic model based on his mathematical theory 
was developed in the field of financial technology and it 
is widely accepted in the real financial market. 

In the 1990s, analysis of high frequency tick-by-tick 
data by physicist clarified that market price is not a sim- 
ple random walk[l3j|. The market price has some em- 
pirically stylized facts which clearly deviate from a pure 
random walk. In order to build a model which fullfills 
those properties, many variants of random walk models 
of market prices have been proposed. 

Recently, one of the authors (M. T.) and her group in- 
troduced a new type of market model called potentials of 
unbalanced complex kinetics (PUCK), wherein the mar- 
ket price is described by a random walker in a tempo- 



rally deforming potential force whose center is given by 
the trace of the walker 0, HE 0- This model has 
quite different statistical properties from the case of the 
fixed potential function and its center, i.e., the Ulenbeck- 
Ornstein process 17|. Moreover, the continuum limit of 
PUCK model is equivalent to the Langevin equation with 
time-dependent coefficients [l8j|. 

From the viewpoint of this model with a quadric po- 
tential function, market states are categorized into five 
conditions, characterized by the potential function: (1) 
a pure random walk state, which is given by the case of 
no potential force. In mathematical finance, this condi- 
tion is assumed to the generic property of the price fluc- 
tuation of a financial market; (2) a stable state, which 
is described by attractive potential force. Peinke and 
his group show existence of similar kind of attractive 
potential in foreign exchange market Ausloos and 
his group applied the same method and obtained simi- 
lar results for the stock index data[2(|; (3) an unstable 
state which is described by a repulsive potential force; 
(4) a non-stationary state, which is characterized by a 
strong attractive force causing the exponential oscilla- 
tion of price plj]; and (5) a non-stationary state, which is 
characterized by a strong repulsive force causing nearly 
monotonic price changes as a result. Data analysis based 
on PUCK has clarified that the market states are chang- 
ing in various timcscalcs. 

The PUCK model is known to satisfy basic market 
statistics such as the power-law distribution of price 
changes [H, [23|, [HI, rapid decay of the autocorrela- 
tion of price changes, long correlation of the square of 
price changes fl3|. and abnormal diffusion properties [2rjj 
in a short timescale 27[. It has been clarified that 
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the market potential force is closely related to the mass 
behavior of dealers, especially, the trend- following be- 
havior [28| . Morever, the Nobel prize laureled market 
model - the autoregressive conditional heteroskedasticity 
(ARCH) model - is derived in a very special limit case, 
where a nonlinear potential force appears with randomly 
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chosen signs at each time step [29[. 

In this paper, we reformulate the data analysis method 
of PUCK taking into account the higher-order potential 
functions, and we focus our attention on finding precur- 
sors of bubbles and crashes. The data we analyze here as 
an example is the tick-by-tick data of yen-dollar exchange 
rates, which showed the largest rate change in 1998. 



II. THE PUCK MODEL 

The PUCK model is represented by the following set 
of equations. 

d 



P(t + 1)-P(t) = - U(p,t)\ p=P{t) _ PM{t) 
+/(*) 

1 Af-l 



(1) 
(2) 



k=0 



Here, P(t) is the noise-reduced market price at the t- 
th tick obtained by applying the optimal moving aver- 
age (30j . f(t) is the random noise typically a Gaussian 
white noise, U (p, t) is the potential function which is ap- 
proximated by the following Taylor expansion form with 
time-dependent coefficients bk(t), 



(3) 



fc=l 



Pm it) is the average of the past price changes with size M 
ticks, and it is assumed that this point gives the center 
of the potential function that moves with the random 
walker's footprint (Fig. [T]). For the yen-dollar market a 
typical value of M is smaller than ten. We postulate that 
it is an averaged expectation price of the market dealers. 

In Fig. [2] we show typical examples of potential func- 
tions observed in the yen-dollar market for given 2000 
data points. Here, the mean interval of tick is about fif- 
teen seconds in both cases. The market fluctuations of 
Fig. an d Fig. [2Ji look very different intuitively; how- 
ever, the square of the market price fluctuations in unit 
time, called the volatility, takes about the same value in 
both cases. By applying the PUCK model we can clearly 
observe the difference in the quadric term of the potential 
function as shown in Fig. [2J; and Fig. [2f. In the case of 
Fig. [2t the value of 62 it) is positive and the market fluc- 
tuation is stable, while in the case of Fig. [2f the value of 
&2(i) is negative and market fluctuation is unstable. As 
typically shown here the potential coefficients other than 
62 (t), such as bi(t) or 63 (t), are generally negligibly small 
in ordinary market states. It is found that these quadric 
potential functions can generally be found in ordinary 
states of financial markets; however, in particular peri- 
ods, when prices move rather monotonically we can ob- 
serve a new type of nonquadric potential functions such 
as a cubic or a quartic function, as shown later. 
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FIG. 1: Schematic diagram of the PUCK model. For 
given market time series we calculate P(t + 1) — P(t) and 
P(t) — Pmit) as shown in the upper diagram. The plots in 
the lower diagram is proportional to the derivative of the po- 
tential function, —-^U(p,t). 



In order to describe such nonlinear potential functions 
correctly, we introduce a higher-order term in the poten- 
tial model as follows. 



U( P ,t) 



7 + 1 



7+1 



(4) 



This higher-order potential model can reproduce all 
kinds of typical price motions in bubbles, crashes and 
hyper-inflations as mentioned in the introduction of this 
paper [31 1 . 

When the market potential function is described by a 
cubic function with a local minimum point, there are 
obviously two states for the random walker as shown 
schematically in Fig. [31 The first state is the trapped 
state that is practically the same as the case of a ran- 
dom walk in a stable quadric potential function, however, 
with a finite possibility the walker goes over the poten- 
tial barrier. The second state is the nontrapped state, 
wherein the walker goes down the potential slope indefi- 
nitely following a double-exponential growth or even the 
finite time divergence in the continuum limit. In partic- 
ular, we can expect a sudden transition of the movement 
of prices from a stable state to a highly unstable state 
even when the market potential function is invariant. 



III. ESTIMATION OF NON-LINEAR 
POTENTIAL FUNCTIONS 

For estimating the nonlinear potential function from 
real data of bubbles or crashes, we need to consider the 
meaning of probability distribution of the independent 
random variable f(t) in Eq. (Q]). From a skeptical view- 
point about the existence of the potential forces, one 
might think that any price movement can be fully de- 
scribed by only the random noise term fit), assuming 
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FIG. 2: Examples of estimated potential functions. For given yen-dollar rates, Fig. [2^ and Fig. [2}i, the relations between P(t) 
(heavy line) and P/vr(i) (thin line) are observed as shown in Fig. [2b and Fig. [2^, respectively. The potential functions, Fig. [2}: 
and Fig. [5f, are estimated by integrating the plots of the lower diagram of Fig. Q] 
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"Nontrapped state" "Trapped state" 

FIG. 3: Schematic diagram of the transition of price changes 
from a stable trapped state to a highly unstable nontrapped 
state. 



that f(t) takes extremely asymmetric values repeatedly 
by chance. This case is theoretically realizable; how- 
ever, the probability of the occurrence of the whole event, 
which is estimated by the products of probability density 
of f(t), becomes negligibly small for the cases of bubbles 
or crashes. This is confirmed systematically in the fol- 
lowing manner. 

For the given symmetric probability distribution of 
/(<), w(f(t)), we can calculate the probability of the oc- 
currence of any time series by assuming that the pro- 
cess is governed by Eq. |T]) with the potential function 
given by Eq. (|4]) , and the parameters of the model are 
constants, namely &i(t) = &x,& 7 (i) = & 7 , in the period 
from t = n to n + N. The probability of the occur- 
rence or the likelihood of a given market price time se- 
ries, {P(n),P(n+l), - ■ ■ ,P(n+N)}, is calculated by the 
following equation. 



n+N 



i(h, 7 ,b 7 ,M)= n w (f(t)) 



(5) 



density of f(t), we apply the Gaussian distribution with 
zero mean. 

In statistics, this type of maximum-likelihood method 
is very popular, and the method of choosing the optimum 
number of parameters has already been established. In 
general, a model with more parameters has a greater ten- 
dency to receive a higher value of likelihood. In order to 
solve an overfitting problem and to evaluate which of the 
two models with different number of parameters is bet- 
ter, we should not use the likelihood itself, but apply the 
information criterion called Akaikc information criterion 
(AIC) that is composed of the logarithm of likelihood 
and the term of pen alty, which depends on the number 
of parameters, k [32l. |33|. 



AIC = -21nZ(&i,7,& 7 ,M) + 2fc 



(0) 



where each value of f(t) is determined from Eq. |T]) with 
Eq. ((31). As for the functional form of the probability 



Using this method, the contour plot of the value of AIC 
for a given time series (period II in Fig. [5]) is plotted in 
Fig. [U Here, the control parameters for determining 
the center of potential force are b\\ 7; 6 7 and M, and 
the value of the moving average size M. The case of 
no potential force corresponds to the parameter values 
b\ = 0, and 6 7 = 0. From these figures we can find a non- 
trivial set of parameters that maximize the likelihood. 

We can introduce, for example, one more nonlinear 
term in Eq. (0J and compare the value of AIC with 
that of one nonlinear term model. As a result, the best 
market price model in view of the AIC value is given 
by the model already introduced in Eq. ^ with the 
lowest integer nonlinear exponent 7 = 2, that is, the 
cubic potential function. We also apply the Bayesian 
information criterion (BIC) and replace the distribution 
of f(t) by a distribution with long tails, but the results 
are almost the same in any case. 

In Fig. [5j we illustrate a part of the yen-dollar ex- 
change rates from September 23rd to October 26th of 
1998, during this period the yen-doller market showed 
historically the largest fluctuations. In period I, the po- 
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FIG. 4: Parameter space of the values of AIC and the optimal 
parameters (®) in the case of Fig. [5} II. bi vs bi (a), M vs bi 
(b) 

tential function is quadratic and stable. In period II, the 
left side of the potential function becomes shorter and 
the best fit cubic function indicates the possibility that 
the exchange rate might crash. In period III, we can 
clearly observe the cubic potential with a very shallow 
well, which indicates that the exchange rate may tend 
to decrease monotonically. In period IV, the main crash 
is over, however, the potential function is unstable and 
symmetric, which implies that the market will be turbu- 
lent, but the fluctuations will be symmetric. 

We confirmed that these results do not change even 
if we replace the distribution function of f(t) from the 
normal distribution to a fat-tailed distribution or to an 
asymmetric distribution. It should be noted that we find 
that the potential function becomes cubic a little before 
the largest crash, which occurs in the middle of the figure. 
This indicates that the cubic potential function can be a 
precursor to the largest event. 

IV. CONCLUSION 

In this paper, we introduced a higher-order potential 
model suitable for the analysis of the market forces ob- 
served in financial bubbles or crashes. In order to show 
the statistical significance of the existence of higher-order 
potential functions, we applied the information criterion. 
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